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In this paper we have considered a model of modified Chaplygin gas and its role in accelerating 
phase of the universe. We have assumed that the equation of state of this modified model is valid 
from the radiation era to ACDM model. We have used recently developed statefinder parameters 
in characterizing different phase of the universe diagrammatically. 
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I. INTRODUCTION 

Recent observations of the luminosity of type la supernovae indicate [1, 2] an accelerated 
expansion of the universe and lead to the search for a new type of matter which violate the 
strong energy condition p + 3p < 0. The matter consent responsible for such a condition 
to be satisfied at a certain stage of evolution of the universe is referred to as dark energy. 
There are different candidates to play the role of the dark energy. The type of dark energy 
represented by a scalar field is often called quintessence. The transition from a universe 
filled with matter to an exponentially expanding universe does not necessarily require 
the presence of the scalar field as the only alternative. In particular one can try another 
alternative by using an exotic type of fluid - the so-called Chaplygin gas which obeys an 
equation of state like p — —B/p [3] where p and p are respectively the pressure and energy 
density and B is a positive constant. Subsequently the above equation was modified to 
the form p = —B/p a with < a < 1. This model gives the cosmological evolution from 
an initial dust like matter to an asymptotic cosmological constant with an epoch that can 
be seen as a mixture of a cosmological constant and a fluid obeying an equation of state 
p = ap. This generalized model has been studied previously [4, 5]. 

In the paper [4] , the flat Friedmann model filled with Chaplygin fluid has been analyzed in 
terms of the recently proposed "statefinder" parameters [6]. In fact trajectories in the {s, r} 
plane corresponding to different cosmological models demonstrate qualitatively different 
behaviour. The statefinder diagnostic along with future SNAP observations may perhaps be 
used to discriminate between different dark energy models. The above statefinder diagnostic 
pair are constructed from the scale factor a(t) and its derivatives upto the third order as 
follows: 



r = — — - and s — — -, =-r- (1) 

aH3 3(<z-i) 

where H and q (= — ff) are the Hubble parameter and the deceleration parameter 
respectively. These parameters are dimensionless and allow us to characterize the properties 
of dark energy in a model independent manner. The parameter r forms the next step in 
the hierarchy of geometrical cosmological parameters after H and q. 

In section II, all the discussions are valid in general for k — 0, ±1, but in section III we have 
specifically considered the simple case of a spatially flat universe (k = 0), which naturally 
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follows from the simplest version of the inflationary scenario and is confirmed by recent 
CMB experiments [7]. In our present work we consider a more general modified Chaplygin 
gas obeying an equation of state [8] 



p = Ap-— with < a < 1 (2) 

This equation of state shows radiation era (when A = 1/3) at one extreme (when the 
scale factor a(t) is vanishingly small) while a ACDM model at the other extreme (when the 
scale factor a(t) is infinitely large). At all stages it shows a mixture. Also in between there 
is also one stage when the pressure vanishes and the matter content is equivalent to a pure 
dust. We have further described this particular cosmological model from the field theoretical 
point of view by introducing a scalar field (f> and a self-interacting potential V(<p) with the 
effective Lagrangian 



In the paper of Gorini et al [4], it has been shown that for the simple flat Friedmann 
model with Chaplygin gas can equivalcntly described in terms of a homogeneous minimally 
coupled scalar field (p. In this case FRW equations for Chaplygin gas fits in Barrow's scheme 
[9]. Following Barrow [10], Kamenshchik et al [3, 11] have obtained homogeneous scalar 
field 4>{t) and a potential V{4>) to describe Chaplygin cosmology. In the section II, we have 
obtained the corresponding expressions for a modified Chaplygin gas as generalization of 
the previous work. 



It has been possible to express V{4>) in terms of the scale factor a(t) for arbitrary values 
of the constant a and have explicitly shown how V((j>) varies as the scale factor interpolates 
between radiation and ACDM stages for different values of the constant A. In the next 
stage evolution of the model universe has been studied in the {r, s} plane for the entire 
physically realistic history of the universe and compared to that for either pure Chaplygin 
gas model or a generalized Chaplygin gas model as given in [5]. 



II. MODIFIED CHAPLYGIN GAS IN FRW MODEL 



The metric of a homogeneous and isotropic universe in FRW model is 



ds 2 = dt 2 - a 2 (t) 



dr 2 



+ r 2 (d8 2 + sin 2 0d(t) 2 ) 



1 — kr 2 

where a(t) is the scale factor and k (= 0, ±1) is the curvature scalar. 
The Einstein field equations are 



(4) 



and 



k 1 

- = oP ( 5 ) 



a 



- n =-\(P + ?>P) ( 6 ) 
a o 

where p and p are energy density and isotropic pressure respectively (choosing 
8ttG = c=1). 
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The energy conservation equation is 

P + 'Ap + p) = Q (7) 
Using equation (2) we have the solution of p as 

i 

1 1 ~ 

(8) 



B C 
+ 



1 + A a 3 ( 1 + A )( 1 + a ) 
where C is an arbitrary integration constant. 

Due to complicated form of the equation of state (2), the scale factor a(t) can not be 
solved from the Einstein field equations (5) and (6) for arbitrary k. However for k = 0, the 
solution of scale factor a(t) has the form 

3(1 + A) B 3(1 + A) \/3 

a-^~ 2 F 1 [x, x, 1 + x, - C{1 + A f ~^} = ^-(1 + A)CH (9) 

where x = 2 (i+ a ) an< ^ ^ s ^ ne hypergeometric function. 
Now for small value of scale factor a(t), we have 

^^(T+A)' ( 10 ) 

which is very large and corresponds to the universe dominated by an equation of state 
V = Ap. 

Also for large value of scale factor a(t), 

i 

which correspond to an otherwise empty universe with a cosmological constant (1^4) 1+ ° ■ 
For accelerating universe q must be negative i.e., a > i.e., 

a 3(l + a)(l+A) > C(l + A)(l + 3A) 

IB 

This expression shows that for small value of scale factor we have decelerating universe 
while for large values of scale factor we have accelerating universe and the transition occurs 



when the scale factor has the expression a = 



C(1+A){1+3A) 
2B 



3(l + o)(l + A) 



Considering now the subleading terms in equation (8) at large values of a, we can obtain 
the following expressions for the energy density and pressure: 

B \^C(1 + A\rh o _, (1+a)(1+A) (13) 



1 + AJ 1 + a V B 

and 

B \^ C fl + A\^ 



l+A 1 + a \ B 



{a+(l + a)A} a -3(i+a)(i+A) (14) 
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Equations (13) and (14) describe the mixture of a cosmological constant equal to 



(1^4) 1+a with a matter whose equation of state is given by 

p = {a+ (1 + a)A}p , 
which for a pure Chaplygin gas reduces to p = ap. 



(15) 



Now we consider this energy density and pressure corresponding to a scalar field <f> having a 
self-interacting potential V(<p). The Lagrangian of the scalar field has been given in equation 
(3). The analogous energy density p^ and pressure p^ for the scalar field are the following: 



P4> = \f + V(4>) = P : 



B C 

i + a + a m+A)(i+ a ) 



(16) 



and 



1 J 2 — • B 



P4 , = ^ 2 -V(<j>)=Ap--=A 
z p 



B 



C 



1 + A a 3(l+A)(l+a) 

Hence for flat universe (i.e., k — 0) we have 



-B 



B 



C 



1 + A a 3(l+A)(l+c) 



17) 
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B 



+ 



C 



1 + A a 3 ( 1+A )( 1+a ) 



B 



B 



+ 



C 



1 + A a 3 ( 1+A )( 1+a ) 



and 



V(4>) = -(l-A) 



B 



+ 



c 



1 + A a 3 ( 1+A )( 1+a ) 



l + a i 

+ 2 B 



B 



+ 



C 



1 + A a 3 ( 1+A )( 1+a ) 



(18) 



(19) 



From equation (18) we have the relation between scale factor a(t) and scalar field <fi as 



VS{l + a){l + A) 
From equation (19) we get 



Sinh 



C{l + A) 



B 



a §(l+a)(l+A) 



(20) 



v ^=V i - a A^a 



^ „ ,_a_ f V3\/TTA(1 + a) 
Cosh^ < — -( 



1 / B 
+ 2 B llTA 



1+% goa ^/^±g(l±fOj 



(21) 



Since 2V{<p) = p — p = (1 — A)p + j^. So when p — > oo i.e., when a — > the situations 
are completely different for A ^ 1 and for A = 1. In the second case — > as p — > oo, 

but in the first case V((/>) — > oo as p — > oo. This is why we get qualitatively two different 
pictures in two cases - as reflected in the figures 2 and 5. In the other limit a — > oo, 



V(4>) — > (iTa) ^ or t ^ ie a b° ve two cases. 



For example, A = |, B = C = a = 1, ^(0) — > as a — > oo 
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Figs. 4-6 shows variation of 4> and V against a and <f> for A = 1 and a = 0.6, 1. 



G 



and A = 1, £? = C = a = 1, F(^>) — > ^ as a — > oo. 

The graphical representation of </> against a and V(0) against a and <f> respectively have 
been shown in figures 1-3 for A = 1/3 and figures 4-6 for A = 1 . From figures 1 , 4 we have 
seen that scalar field <j> sharply falls when scale factor a(t) increases, both for A = 1/3 and 
A = 1. In figure 2, we see that potential function V(</>) sharply decreases from extremely 
large value and thereafter it increases to a fixed value for A = 1/3. On the other hand for 
A = 1, the nature of V(4>) against a is as stated in the following. For small a, V(<j>) varies in- 
significantly from zero value and then increases sharply to a constant value. In figure 3, V((f) 
decreases slightly and then increases to infinitely large value for A = 1/3 while for A = 1, 
in figure 6, V(4>) decreases to zero asymptotically starting from a extremely large value. So 
V{<p) demonstrates completely different behaviour for two different values of A = 1 and 
A = 1/3. Further the figures show how V{4>) varies with <f> with different values chosen for a. 



III. THE ROLE OF STATEFINDER PARAMETERS IN FRW UNIVERSE 

The statefinder parameters have been defined in equation (1). Trajectories in the {r, s} 
plane corresponding to different cosmological models, for example ACDM model diagrams 
correspond to the fixed point s = 0, r = 1. 

For Fricdmann model with flat universe (i.e., k = 0), 



and 



a 1 3 p 

So from equation (1) we get 



(23) 



r = l + ?fl + ?^, -=fl + ?^ (24) 



(25) 



2 \ p J dp \ p J dp 

Thus we get the ratio between p and p: 

p _ 2(r- 1) 
p ~ 9s 

For modified Chaplygin gas, velocity of sound can be written as 

From equations (24) and (25) we get the relation between r and s: 

18(r - l)s 2 + 18as(r - 1) + Aa{r - if = 9sA(l + a)(2r + 9s - 2) (27) 

In the {r, s} plane the above equation has only one asymptote parallel to s-axis 
namely, r = 1 + §A(1 + a) and the asymptote intersects the curve in only one point 



Fig.7 



Fig. 7 shows the variation of s against r for different values of a (= 0.6, 1) and for A — 1/3. 



Figure 7 shows the variation of s with the variation of r for ^4=1/3 and for a = 0.6, 1. 
The portion of the curve on the positive side of s which is physically admissible is only 
the values of r greater than |l + |A(l + j4)}. The part of the curve between r = 1 and 
r = 1 + 1 A(l + A) with positive value of s is not admissible (we have not shown that part in 
the figure 7) because for the Chaplygin gas under consideration we face a situation, where 
the magnitude of the constant B becomes negative. It can be shown in the following way. 

Since we consider the range of A as < A < 1/2, we immediately get the inequality 
relation 



1 - 9 A < 



(28) 



for all points lying in the part of the curve between r = 1 and r = 1 + |A(1 + A) as stated 
above. In this range positive values of the parameter s is explicitly given by 



l(r-l)(l-2A) 
2 9A+{l-r) 




{9A+1- r) 
9 (1-2A) 2 



(29) 



The above expression for s is obtained by solving equation (27). One must note that the 
other root containing negative sign before the square root symbol is to be left out since we 
consider only the positive magnitude of s. 



Again from (25) and (2), we get 



V - = A- 



B 2(r - 1) 



1 a+l 



9s 



which in view of (29) leads to the relation 



B 

nO+l 



1 (1 - 2A) 



8 (9A + 1 - r) 

9 (1-2A) 2 



(30) 
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One can show by a simple algebraic exercise that the right hand side of (30) is clearly 
negative for all points corresponding to r < 1 + §A(1 + A) and hence B becomes negative, 
which is unacceptable in view of the Chaplygin gas model we are considering. This upper 
limit for r corresponds to r = 3 for A = 1/3 that is when we choose our universe to be a 
mixture of radiation and Chaplygin gas. 

Thus the curve in the positive side of s starts from radiation era and goes asymptotically 
to the dust model. But the portion in the negative side of s represents the evolution from 
dust state (s = — oo) to the ACDM model (s = 0). Thus the total curve represents the 
evolution of the universe starting from the radiation era to the ACDM model. 



IV. DISCUSSIONS 

In this work, we have presented a model for modified Chaplygin gas. In this model we 
are able to describe the universe from the radiation era (A = 1/3 and p is very large) to 
ACDM model (p is small constant). So compare to Chaplygin gas model, the present model 
describe universe to a large extent. Also if we put A = with a = 1, then we can recover 
the results of the Chaplygin gas model. In figure 7, for {r, s} diagram the portion of the 
curve for s > between r — 1 tor = 1 + §A(1 + A) is not describable by the modified 
Chaplygin gas under consideration. For example, if we choose r — 1.03, A — l/Z then from 
the curve s = 0.01 which corresponds to q = 3/2 and hence we have from the equation of 
state, B < which is not valid for the specific Chaplygin gas model considered here. At 
the large value of the scale factor we must have some stage where the pressure becomes 
negative and hence B has to be chosen positive. It follows therefore that a portion of the 
curve as mentioned above should not remain valid. 
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